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Q '■ ABSTRACT 

V- j ■ X-ray surface brightness fluctuations in the core (650 x 650 kpc) region of the Coma 

c/5 . cluster observed with XMM-Newton and Chandra are analyzed using a 2D power 

^ ■ spectrum approach. The resulting 2D spectra are converted to 3D power spectra of 

gas density fluctuations. Our independent analyses of the XMM-Newton and Chandra 
. observations are in excellent agreement and provide the most sensitive measurements 

^ ■ of surface brightness and density fluctuations for a hot cluster. We find that the char- 

V/^ I acteristic amplitude of the volume filling density fluctuations relative to the smooth 

■ underlying density distribution varies from 7-10% on scales of ^^500 kpc down to -^5% 
OO ■ at scales ~ 30 kpc. On smaller spatial scales, projection effects smear the density fluc- 

, tuations by a large factor, precluding strong limits on the fluctuations in 3D. On the 

■ largest scales probed (hundreds of kpc) , the dominant contributions to the observed 

fluctuations most likely arise from perturbations of the gravitational potential by the 
two most massive galaxies in Coma, NGC4874 and NGC4889, and the low entropy 
gas brought to the cluster by an infalling group. Other plausible sources of X-ray 
surface brightness fluctuations are discussed, including turbulence, metal abundance 
variations, and unresolved sources. Despite a variety of possible origins for density 

■ fluctuations, the gas in the Coma cluster core is remarkably homogeneous on scales 
H ! from - 500 to ~30 kpc. 
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1 INTRODUCTION 

The hot intracluster medium (ICM), which fills the cluster 
gravitational potential, is the dominant baryonic component 
of rich galaxy clusters. In a smooth and static gravitational 
potential, the gas in hydrostatic equilibrium should have a 
smooth density and temperature distribution aligned with 
the equipotential surfaces. From X-ray data both the density 
and temperature radial profiles are routinely measured, and 
the comparison of mass estimates from hydrostatic equilib- 
rium and lensing suggests overall good agreement (to within 



10-20%). However, there are a variety of reasons why the gas 
properties may be perturbed on a range of spatial scales: 
perturbation of the gravitational potential, turbulent gas 
motions, imperfect mixing of the gas with different entropies 
displaced by gas motions, presence of bubbles filled with rel- 
ativistic plasma, etc. These perturbations of the ICM prop- 
erties can be detected in raw observational data as flux (sur- 
face brightness) fluctuations, variations of temperature, or 
variations of the projected velocity field (which can be mea- 
sured with future high resolution X-ray calorimeters, like 
ASTRO-H). 
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Several previous studies have explored fluctuations in 
galaxy clusters. X-ray surface brightness fluctuations on ar- 
cminute scales, relative to a smooth model, were detected 
for a sample of clusters us ing the Einstein Observatory data 
(|Soltan fc Fabricant|[l990l '). Projected pressure fluctuations 
in the Com a cluster, using XMM-N ewton observations were 
studied bv ISchuecker et alj (12004 ). Another approach was 
adopted by Kawahara et al.l ( 20081 ') . where the distribution 
function of the surface brightness fluctuation was calculated 
for the A3667 cluster and linked to the 3D density fluctua- 
tions. 

Here we present an analysis of the X-ray surface bright- 
ness fluctuations in the core of the Coma cluster using 
both Chandra and XMM-Newton observations, calculate the 
power spectrum of these fluctuations and estimate the am- 
plitudes of the 3D density fluctuations on scales ranging 
from 30 kpc to 500 kpc. T he projected pre s sure fluctua- 
tions have been analyzed by ISchuecker et al.l (|2004l l over a 
limited range of spatial scales. Using X-ray surface bright- 
ness instead of pressure allows one to probe smaller spatial 
scales and greatly simplifies accounting for statistical noise. 

In our work the choice of the Coma cluster was moti- 
vated by three factors: proximity (thus allowing us to probe 
small physical scales) , presence of a large flat surface bright- 
ness X-ray core (simplifying the conversion of 2D power 
spectra to 3D) and brightness of the cluster. The exten- 
sion of the analysis to even brighter and closer clusters, in 
particular, Perseus and M87, is straightforward and will be 
presented elsewhere. 

Throughout the paper the redshift of Coma was as- 
sumed to be z = 0.023, corresponding to an angular di- 
ameter distance of 93 Mpc (for h = 0.72); 1" corresponds to 
0.45 kpc. 

This paper is organized as follows: in ij2]we describe the 
data-sets used in the analysis. In i|3]the relation between 2D 
and 3D power spectra is explained. The measurements of 
the density fluctuation power spectrum are presented in ^ 
The implications of our results are discussed in JS] and 311 
contains our conclusions. 



2 DATA AND INITIAL PROCESSING 

For our analysis we used publicly available XMM-Newton 
and Chandra data. 

For Chandra, the following OBSIDs were used: 555, 556, 
1086, 1112, 1113, 1114, 9714, 10672. The data were p r epare d 
following the procedure described in lVikhlinin et aU (|2005h . 
This includes filtering of high background periods, appli- 
cation of the latest calibration corrections to the detected 
X-ray photons, and determination of the background inten- 
sity in each observation. The total exposure of the resulting 
data set is ~ 1.15 x 10^ s. 

For XMM-Newton, we used a large set of pointings, 
covering a field of more than 1x1 degrees (although only 
the most central part of this field was used in the subse- 
quent analysis). The data were prepared by removing back- 
ground fiares using the light curve of the detected events 
above 10 keV and re-normalizing the "blank fields" back- 
ground to match the observed count rate in the 11-12 keV 




Figure 1. XMM mosaic image of Coma (1° X 1°, or 1.6 X 1.6 
Mpc). Boxes show the regions used for the power spectrum cal- 
culation: central 14' X 14' box - Chandra field; larger 24' X 24' 
box - XMM-Ncwton field. The smallest 13' X 13' box near the top 
of the mosaic, which contains the X-ray bright Seyfert galaxy X- 
Comae, was used to determine the shape of the power spectrum 
produced by point sources in XMM-Ncwton images. 



band. In the subsequent analysis we used the data from the 
EPIC/MOS (European Photon Imaging Camera/Metal Ox- 
ide Semi-conductor) detector. The total exposure of all used 
observations is 880 ksec, covering a much larger area than 
the Chandra observations (not only the cluster core but also 
the outskirts of the cluster). 

The images in the 0.5-2.5 and 0.5-4 keV bands were gen- 
erated for XMM-Newton and Chandra respectively, along 
with the exposure and background maps, which were used 
for subsequent analysis. 

The background-subtracted, exposure-corrected XMM- 
Newton image (1° x 1°, or 1.6 x 1.6 Mpc) is shown in 
Fig[T] For the subsequent analysis of XMM-Newton data 
only the central 24' x 24' part of the image was used (roughly 
corresponding to the flat surface brightness cluster core). 
For Chandra we also excluded edges of the image with low 
exposure, keeping only the central 14' x 14' region. 

An azimuthly averaged X-ray surface brightness pro- 
flle I{R) is shown in Fig[2] The observed distribution 
can be reasonably well describe d by a simple /3-model 
(|Cavaliere fc Fusco-Femiandll978l ) with core radius Tc = 9' 
(or ~245 kpc) and (3 = 0.6 

m = ^°„. 3/^-0.5 . (1) 



1+ f 



where R is the projected distance from the cluster center 
and 7o is the surface brightness at the center. In i]4.5l we 
use this model (and a more complicated one) to remove the 
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Figure 2. Radial X-ray surface brightness profile (XMM- 
Newton) of the Coma cluster. The solid line shows the surface 
brightness for a /3-model with core radius rc = 9' (or ~245 kpc) 
and 13 = 0.6. 



global cluster emission. As discussed in ^5. li the results of our 
analysis are only weakly sensitive to the particular choice of 
the /3-model. Partly this is because our analysis is apphed 
to the core region of the Coma cluster. 



3 3D TO 2D PROJECTION 

We assume that the density distribution in the cluster is 
described by a simple decomposition into undisturbed and 
fluctuating components 



ni^, y, z) = no{x, y, z) [1 + 5[x, y, z)] , 



(2) 



where no{x, y, z) is the undisturbed density distribution (e.g. 
standard spherically symmetric /J-modeJj shown in Figia 

and the dimensionless quantity 5{x,y,z) = — ^ 1 de- 
no 

scribes the fluctuating part. For ICM temperatures above 
~3 keV, the emissivity of plasma in the soft band, which 
dominates the flux detected by XMM-Newton and Chandra, 
mainly depends on the square of the electron density and is 
only weakly sensitive to the temperaturfl Accordingly the 



^ One can imagine more complicated models, e.g. ellipsoidal j3- 
model or a model with more sophisticated radial density profiles. 
In J5] we argue that the freedom in choosing the undisturbed 
model mostly affects the power spectra at large scale. 

The impact of possible abundance variations is briefly discussed 



volume emissivity e can be written as 

£{x,y,z) = Cn'^(x,y,z) Cnl{x,y,z) [I + 2S{x,y,z)\ , (3) 

neglecting terms of order 5^ and C is the coefficient relat- 
ing the square of the density and the emissivity. Since we 
are interested in fluctuations, relative to the average sur- 
face brightness profile, for the subsequent analysis we set 
C = 1. We further assume that 5{x,y,z) is a homogeneous 
and isotropic random field, characterized by the power spec- 
trum P-iD- 



Psoikx, ky, k^) = 



S{x,y,z)e 



dxdydz 



(4) 



By isotropy assumption, PiD{kx,ky,kz) — P3D(\k\), where 

The observed image I{x, y) is the integral of the emis- 
sivity along the line of sight: 



/ 



I{^,y) = / e{x,y,z)dz = 



(5) 



^(a;,y) + 2 / no{x,y,z)5{x,y,z)dz, 



where Io{x,y) is the undisturbed surface brightness. Divid- 
ing the observed image by Io{x,y) yields the normalized 
image J{x, y) 



J{x,y) = 1-^2 



S nl{x,y,z)5{x,y,z)dz 



J nl{x,y,z)dz 
The 2D power spectrum of J{x, y) 



P2D{kx,ky) = 



J{x,y)e 



-i27r{xkx-\-yky ) 



dxdy 



(6) 



(7) 



can be obtained directly from observational data. By 



isotropy, P2D = P2D{kxy), where k^y = ^kl +k^. The ex- 
pression ([6]) contains the following function of x, y and z 
under the integral 

no{^,y,z) _ nl{x,y,z) 



J nl{x,y,z')dz' Io{x,y) 



(8) 



For the range of projected distances we consider, the weak 
dependence of rj{x, y, z) on x and y can be neglected (the 
changes in ri{x,y, z) are at most factor of 2 across the whole 
image 24' x 24' for the adopted /3-model). Thus r](x,y,z) ~ 
r][z) and the P2D spectrum can be link ed to the P^p of 
5{x, y , z) as, follows (|Peacockll 19991 . see also lZhuravleva et al.l 
I2OI1I ) 



{k^y) = ^ PsD [^Jkly + ki^ \W{k,)\'dk. 



(9) 



where \ W{k^)\'^ \s the ID power spectrum of the normalized 
emissivity distribution along the line of sight. 



W{k, 



r][z)e ~az. 



(10) 



Note that the properties of the emission measure dis- 
tribution along the line of sight determine the characteristic 
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wavenumber kz,cutoff above which |M^(fc2)|^ falls off towards 
larger k. For k kz, cutoff the expression ([9]) simplifies to 



P2Dik) ^ 4P3Dik) J \W{k.)fdk, 



(11) 



i.e., the 2D power spectrum of the surface brightness fiuctu- 
ations is essentially equal to the 3D power spectrum of the 
density fluctuations apart from the normalization constant 
4 J \ W{kz)\'^dkz, which is easily measured for a cluster. In 
principle J \ W{ki;)\'^dkz depends on x and y and so is differ- 
ent for different lines of sight, but as we show in FigO for 
the central region of Coma this difference can be neglected. 

For practical reasons it is easier to use the character- 
istic amplitude of the fluctuations, rather than the power 
spectrum. The amplitudes corresponding to the 3D and 2D 
spectra are deflned as 



Asoik) = y^PiD{k)4Tvk^, 



Aioik) = y/P2D(k)2Tvk^ 



(12) 



(13) 



To illustrate the above relations for the case of Coma, 
we set the density distribution eq.Q to the /3-model with 
P = 0.6 and core radius rc — 245 kpc corresponding to 
the azimuthally averaged Coma surface brightness profile, 
shown in Fig[2] 



no(r) oc ( 1 + 



(14) 



and calculated the power spectra of the line-of-sight emis- 
sivity distribution at several projected radii. We further 
parametrize the 3D power spectrum of density fluctuations 
with a simple cored power law 



Psoik) = 



B 



1 + 



fc2 



1/2 ' 



(15) 



where a is the slope of the power spectrum at large k, where 
km is the characteristic wavenumber below which the 3D 
spectrum is fiat and B is the normalization constant. This 
choice of the power spectrum parametrization is rather ar- 
bitrary. It is used only to compare the exact calculation of 
P2D{k) using eq.® with the approximate formula eq. (|lip . 

Using eq. (|15p and the power spectrum |W(fcz)|^, the 
power spectrum of the surface brightness fiuctuations can 
be easily calculated via eq. ([9]). The resulting characteristic 
amplitudes Asn and A2D for different values of km, a and 
different lines of sight (projected distances) are shown in 
Figs. [3l [4] and [5] respectively. 

In each of these images the dotted lines show the ra- 
tio of the 2D and 3D amplitudes, calculated using the exact 
relation ((9| divided by the ratio evaluated using the sim- 
plified relation (|ll|l . which corresponds to the case of small 
scale perturbations. Clearly, independently of the shape of 
the power spectrum (Fig[3]and Fig|4]) or projected distance 
(Fig[5]) the simplified relation (|lip provides an excellent ap- 
proximation of the relation between 2D and 3D spectra (un- 
less the power law spectrum extends to very low k, corre- 
sponding to scales of a few Mpc). Furthermore for a given 
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Figure 3. 3D and 2D amplitudes as a function of k for a model 
power spectrum and for a /3-model cluster with /3 = 0.6 and 
rc = 245 kpc. A cored power law model is used for the 3D 
power spectrum, according to eq.|[l5}. Different colors correspond 
to different break wavenumber (in kpc~^) in the 3D power spec- 
trum: km = 3 W-^ (blue), 3 lO"" (green), 3 10"^ (red), 3 IQ-^ 
(black). Projection (see eq.jgj) is done along the line-of-sight go- 
ing through the Coma center. Solid curves correspond to 
dashed - to A2D, dotted - to the ratio of A2D to the prediction 
of the simplified expression Ullj l. All dotted curves are very close 
to unity at fc ^ 10~^ kpc~^, indicating that expression lllll l can 
be used to relate 3D and 2D amplitudes for k ^ 10~^ kpc~^. In 
this plot, the factor 4 appearing in eq. (|9j is omitted. 



3D spectrum, the expected 2D spectra, measured at differ- 
ent projected distances from the Coma cluster, differ only in 
amplitude, and this difference is less than a factor of 2 (pro- 
vided that the projected distance is less than ~ 600 kpc). 



4 MEASURING THE 3D POWER SPECTRUM 
OF DENSITY FLUCTUATIONS 

Our goal is to obtain 3D power spectra of the density fluctu- 
ations. Our general strategy is as follows. We start with the 
analysis of the 2D power spectra of raw images and evalu- 
ate the contribution of the Poisson noise and point sources. 
Then we pick a smooth model as a representation of the 
undisturbed cluster density distribution, predict the corre- 
sponding surface brightness distribution and calculate 2D 
fluctuations relative to this model. The resulting 2D power 
spectra are corrected for the telescope Point Spread Func- 
tion (PSF) effects. We then convert 2D power spectra to 
3D spectra using eq. (|ll|l and recast the results in terms of 
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Figure 4. 3D (solid) and 2D (daslied) amplitudes as a function 
of k for a model power spectrum. Different slopes of the 3D power 
spectrum a =3.5, 4, 4.5 are shown with blue, red and black col- 
ors respectively, km = 3 10"'' kpc"'^ is the same for all plots. 
The dotted line shows the ratio of A2D to the prediction of the 
simplified expression 1111^ . 



Figure 5. 3D (solid) and 2D (dashed) amplitudes as a function 
of k for a model power spectrum. The 3D power spectrum has 
km = 3 10~^ kpc~^ and a = 4. The 2D amplitude is calculated 
for different projected distances _R =0, 300, 435, 600 kpc (red, 
blue, green, black respectively). The dotted line shows the ratio 
of A2D to the prediction of the simplified expression 



3D amplitudes of density fluctuations relative to the smooth 
underlying 3D model. 



4.1 2D power spectrum evaluation 

The method used to evaluate low resolution power spec- 
tra of the surface brightness fluctuations is described in 
Arevalo et al.l I2OIII (see also lOssenkopf. Krips. fc Stutzkil 
2008h . Briefly one convolves the image with a set of 
Mexican-Hat filters ( with different spatial scales) F{x) — 



X ^ n 

1 5- e~ 2^ , where a is the characteristic spatial scalqj, 

(7^ J 

and relates the variance of the resulting image to the power 
at a given scale. The presence of a non-uniform exposure 
map or gaps in the data (e.g. excised regions around bright 
point sources) is accounted for by representing the Mexican 
Hat filter as a difference between two Gaussians with only 
slightly different width. For each sp atial scale g, the foU ow- 
ing simple steps are performed (see lArevalo et al]l201ll for 
a detailed description): 

(i) the observed image I and the exposure map E axe 
smoothed with two Gaussians G having different smooth- 



^ Mexican-Hat filter has a positive core and negative wings, and 
preferentially selects perturbations with a given scale. 



ing lengths cri = a j \/\ + e and (T2 = (t\/\ + e, where e <C 1 
is a dimensionless number; 

(ii) smoothed images are divided by smoothed exposure 
maps and regions with low (unsmoothed) exposure are 
masked out; 

(iii) the difference of the two images is calculated, this dif- 
ference image is dominated by fluctuations at scales ~ a; 

(iv) the variance of the resulting image is calculated and 
re-cast into an estimate of the power, using a normalization 
coefficient that depends on e. 



Thus the variance is calculated for an image 



I„ = Ex 



G(ji ^ iraw G<y2 ^ ^rau 



(16) 



where E — E{x,y) is the exposure map, Eam = I{x,y) is 
the raw image in counts and Ga o / denotes a convolution 
of an image I with a Gaussian. The final expression for the 
power at a given spatial scale (or given wavenumber) is then 



P2D{k) 



(17) 



wher e summation i s ove r the image pixels (see appendix 
1 m lArevalo et al] [ioij). The resulting 2D power spec- 
trum for the XMM-Newton image of Coma is plotted in 
Fig[6] (dashed line). Effectively the power spectrum shown 
in Fig|6] is a convolution of the true underlying power spec- 
trum with a broad {Ak ~ k) filter in k space. Therefore 
individual measurements of the power are not independent, 
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Figure 6. Raw power spectrum P2D{k) (dashed line), eq.l ll7l l. 
of the XMM-Newton field with sources and Poisson noise. The 
Poisson noise level, estimated using eg.l llSII . is shown with the 
thin horizontal line. The Poisson noise evaluated by adding extra 
Poisson noise to the images (see i|4.2ll is shown with circles. The 
Poisson-subtracted power spectrum is shown with the solid line 
and black squares. The errors shown are due to Poisson noise only 
and do not account for the stochastic nature of the signal itself. 



unlike components of the conventional Fourier transform 
in a rectangular region. However the procedure adopted 
here has a clear advantage compared to the conventional 
Fourier transform, since most of the problems caused by 
non-periodicity of the data-set and the presence of the gaps 
in data are automatically correcte d for (|Arevalo et al.|[201ll : 
lOssenkopf. Krips. fc StutekHbOOSh . 

For XMM-Newton and Chandra observations, we de- 
cided to use images without instrumental and cosmic back- 
ground subtraction. Inclusion of the background does not 
significantly affect the results, but in the case of XMM- 
Newton it would add some high frequency noise due to weak 
compact sources, which were not eliminated completely in 
the blank-field data-sets. 



4.2 Poisson noise 

The Poisson noise contribution can be easily estimated for 
the raw images. We expect a white noise *(flat) power spec- 
trum* for this component with the power: 

Poisson — Si?*^ ' V / 

where summation is over the image pixels. The result of this 
calculation is shown in Fig[6]with the solid horizontal line. 



Very small deviations from this simple estimate are pos- 
sible due to non-trivial behavior of the filter at the edges of 
the image or near the data gaps (partly due to the coarse 
representation of the Mexican Hat filter on a grid for the 
highest frequency bins). To take care of these possible devi- 
ations we made 100 realizations of fake raw images by adding 
extra Poisson noise to each pixel with the mean value set by 
the actual number of counts in the pixel. The whole pro- 
cedure of evaluating the power spectrum is then repeated 
and the mean of the resulting power spectra is calculated. 
Effectively this mean power spectrum has an extra Poisson 
noise component (compared to the original image) with all 
properties similar to the Poisson noise in the original image. 
Therefore the difference between the mean power spectrum 
and the original power spectrum provides an estimate of the 
Poisson noise contribution, which includes all possible effects 
of the filter implementation. 

As expected, the derived noise spectrum has an essen- 
tially flat spectrum. The same procedure allows us to evalu- 
ate uncertainties caused by the presence of Poisson noise in 
the power spectrum. The estimated level of the Poisson noise 
and corresponding uncertainties (in evaluating the Poisson 
noise) are shown in Fig[6]as black circles with error bars. 

The power spectrum with the Poisson noise contribu- 
tion subtracted is shown with the solid curve. The errorbars 
shown are due to the Poisson noise only. 



4.3 Point sources 

The Poisson-noise-subtracted power spectrum (solid line in 
FigO changes slope at fc ~ 0.01 arcsec"^. At larger k, the 
spectrum flattens. This flattening is at least partly caused 
by the presence of compact sources in the image. Indeed, ig- 
noring the blurring of an image by the PSF of the telescope, 
the presence of compact sources (i.e. fsrcS{x — Xsrc,y — ysrc), 
where 5 is the Dirac delta-function at the position of a 
source with flux fsrc), randomly scattered across the im- 
age, produces a flat power spectrum over the entire range 
of wavenumbers. Accounting for the telescope PSF, having 
Gaussian width a, the expected contribution of point sources 
has the shape of the PSF power spectrum, which is flat at 
k <^ 1/a and falls off at larger wavenumbers. 

As a first step towards removing the contribution of 
compact sources we excised a 15" circle (radius) from the im- 
ages around each bright, individually detected source. The 
exposure map has been modified accordingly. The power 
spectra calculated from the images without bright sources 
are shown in Fig[7] with the blue curves. Clearly, much of 
the high-frequency part of the power spectrum goes away 
once bright sources are excluded. For comparison, the raw 
spectrum is shown with black symbols. 

To confirm still further that much of the detected 
higher- frequency power is attributable to the point sources, 
we computed the power spectrum for images containing a 
large number of bright compact sources. 

For Chandra, we used one of the long observations of 
the Lockman Hole with ACIS-I and used a 14' x 14' patch 
of the image. The power spectrum obtained from this field 
after correction for Poisson noise is shown in FigO We also 
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Figure 7. Power spectra of raw images with point sources (black solid) and with 15" circles around point sources excluded (blue solid) 
for XMM-Ncwton field (left) and Chandra field (right). Black dashed lines show the expected shape of a power spectrum associated with 
point sources, taken from the X-Comae and Lockman fields for XMM-Newton and Chandra respectively (see i|4.3l l. Dotted red lines show 
the effect of subtracting an extra component, associated with weak point sources, which are too faint to be detected individually in the 
images. The normalization of this subtracted component is estimated to be ai3% of the power spectrum due to bright detected sources 
(see ^43t . 



plot the same power spectrum in Fig[7]with a dashed line, 
changing only the normalization of the spectrum. Clearly 
the shape of the shoulder in the power spectrum calculated 
for a raw image of Coma (including compact sources) is con- 
sistent with the power spectrum observed in the Lockman 
Hole data. 

For XMM-Newton, we used a patch of the image to 
the NE from the Coma center (Fig[l| . This patch contains 
the X-ray bright Seyfert galaxy X- Comae and a number of 
weaker sources. The power spectrum calculated for this im- 
age is shown in Fig[9l Similarly to Chandra, this power spec- 
trum has a shape consistent with the shape of the shoulder 
obtained from the raw image of Coma (including compact 
sources) as shown in FigO The properties of the PSF vary 
across the telescope field of view. We used a field close to 
the Galactic Plane (OBSID: 0605580901), containing many 
compact sources across the field, to verify that the shape of 
the PSF power spectrum obtained from the X-Comae field 
is not too strongly dominated by the X-Comae itself, which 
was observed close to the center of the XMM-Newton field 
of view. The resulting power spectrum of the Galactic Plane 
image is consistent with the X-Comae patch over the inter- 
esting range of wavenumbers. 

An additional question arises in connection with XMM- 
Newton data; is a 15" circle around a source large enough to 
sufficiently suppress the contribution to the Power Spectrum 
of the PSF wings? To address this question, we repeated 



the calculation of the power spectrum for the patch near 
X-Comae, excluding 15" and 30" circles around each source 
(see Fig llOf) . The corresponding power spectra are shown 
in Fig[^ From the comparison of these spectra with the 
full power spectrum of the X-Comae patch, it is clear that 
removing 15" circles suppresses the high frequency noise by 
a factor of order 300, while 30" circles suppress the signajj 
by a factor of a few 1000. For our purposes, the suppression 
factor of more than a 100 is enough and we fixed the size of 
the excised region around each point source to 15" (radius). 

Apart from the bright detected sources there must be 
faint compact objects that we failed to detect. Their collec- 
tive contribution can be estimated assuming that we know 
the LogN-LogS and limiting fluxes. To this end we calcu- 
lated fluxes of our detected sources using a 15" circle as a 
source region and a twice larger circle as a background re- 
gion for each source (for Chandra the size of the region is 
probably larger than necessary, but it is not very important 
for the present study). The LogN-LogS curves for both in- 
struments are shown in Fig llll (note, that the XMM held is 
larger than the Chandra field). The curves can be reason- 
ably well approximated by the law N oc . Assuming that 



* Note that there is significant contribution of Coma cluster dif- 
fuse emission to the field. For this reason the signal does not go 
to zero even for the larger size of the excised regions. 
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Figure 10. Patch of the XMM-Newton image (13' X 13') near the bright Seyfert galaxy X Comae. In the middle and right panels the 
circles with radius 15" and 30" are excised around bright sources. 
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provided F2 3> i^i, where F2 is the flux of the brightest 
detected source (assuming that the same slope of LogN- 
LogS is applicable to the entire range from Fi to ^2)- For 
XMM-Newton F2 ~ 3000 counts. Thus the contribution 

of unresolved sources can be estimated to be at the level 

Fi 

of — ~ 0.03 relative to the known contribution of bright 
F2 

sources to the power spectrum. We used the difference be- 
tween the power spectra with and without sources to cal- 
culate the power spectrum associated with bright sources 

Pbright — Pwithsrc — Pnosrc aud multiphcd Pbright by 0.03. 

The subtraction of this extra component modifies the high 
frequency side of the power spectrum (see FigO, but the 
magnitude of the modification does not exceed the uncer- 
tainties in the power spectrum, associated with pure Poisson 
noise. Given the uncertainties in estimating this component 
and taking into account that it is mostly important in the 
range of wavenumbers where Poisson noise dominates the 
signal, in the analysis below we neglect this component. 



0.001 0.01 0.1 

k, arcsec ' 



Figure 8. Power spectrum (black) of the Chandra Lockman Hole 
image. The red curve shows the crude approximation of the power 
spectrum of the Lockman Hole field by a simple model (see i|4.4l 
and eql22j 

the same slope continues at lower fluxes we can estimate the 
contribution of faint point sources to the power spectrum: 

Pfa^nt(X / F^ dF (X Fl , (19) 

Jo al< 

where F\ is the minimal flux of the detected source. For 
XMM F\ ~ 100 counts. It is convenient to relate Pfaint to 
the combined contribution of bright detected sources: 

/^'2 fj AT 
—F''dF<xF2, (20) 



4.4 Correction of the Coma power spectrum for 
PSF 

Analysis of the images containing bright point sources can 
also be used to evaluate the effective power spectrum of the 
telescope PSF. Indeed, the convolution of the true sky image 
with the telescope PSF is equivalent to the multiplication of 
the true sky power spectrum by the PSF power spectrum. 
Note that the PSF shape may vary across the instrument 
Field of View (FoV), so one needs to perform the analysis 
separately for different parts of the FoV, or use some ef- 
fective PSF (e.g. FoV averaged). We decided to follow the 
second approach and simply used the observed power spec- 
tra for the Lockman Hole (Chandra) and X-Comae region 
(XMM-Newton) as the PSF power spectrum (see Fig |8l9)l 
approximated with the following simple analytic models: 

PpSF.Chandra = 7 — fT (*^-'-) 

[i+(o^)T 
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Figure 9. Power spectra for a patch near X Comae. Black curve - 
the power spectrum from the image with sources (sec Fig llOl left 
panel). Blue - power spectrum when 15" circles (radius) around 
bright sources are excluded (see Fig llOl middle panel). Magenta 
- power spectrum when 30" circles excluded (see Fig llOl right 
panel) . Thick red curve shows a simple approximation of the PSF 
power spectrum (see i|4.4l and eQ l22l l . Lower red curves show the 
same spectrum scaled down by a factor of 100 and 1000. From 
this plot one can conclude that excising 15" circle around a bright 
source reduces the contribution of the source to the power spec- 
trum by a factor of ~300 or more (see i|4.3l l. 



PPSF,XMM = — 176, (22) 

[l + (o^) ] 

for Chandra and XMM respectively, where the wavenumber 
k is in units of arcsec"^. The correction for PSF blurring is 
done by dividing the power spectrum of the raw image by 
the power spectrum of the PSF. 



Figure 11. LogN-LogS for detected sources in XMM and Chan- 
dra fields (upper and lower red histograms, respectively). The flux 
(in counts) for each source was calculated in 15" circle, using 30" 
for the background. Black solid lines show the law. Assuming 
that the same law continues to lower fluxes, one can estimate the 
contribution of faint (undetected) sources to the power spectrum 
(see ^43t . 

/3-models yield similar power spectra, except at the lowest 
wavenumbers (see i ]5.1|l . 

The resulting 2D power spectra (normalized by the sym- 
metric /3-model, with point sources excised, Poisson noise 
level subtracted) are shown in Fig llSI In this plot the 2D 
power spectrum was converted into the characteristic am- 
plitude of the surface brightness fluctuations relative to the 
smooth model. The measured amplitude of perturbations 
of the surface brightness (before the PSF correction) varies 
from ~10% at scales of ~500 kpc down to less than 1% at 
scales of ~30 kpc. 



4.5 Removing the global surface brightness profile 

As discussed in [JH] it is convenient to divide the surface 
brightness of the raw image by the predictions of a sim- 
ple model, describing the global properties of Coma. The 
divided image can then be interpreted as fluctuations of the 
surface brightness relative to the undisturbed image. To this 
end we experiment with using a simple symmetric /?-model 
(see 321 as an approximation of the cluster global surface 
brightness distribution or a slightly more sophisticated two- 
dimensional /3-model, which allows for two different core 
radii and rotation. For illustration the image divided by 
the symmetric /3-model is shown in the bottom-right corner 
of Fig ll2l In practice both symmetric and two-dimensional 



4.6 Converting to 3D and the amplitude of 
density fluctuations 

The final step of the analysis is the conversion of the 2D 
surface brightness fluctuations (relative to the underlying 
undisturbed image) to the 3D density fluctuation power 
spectrum or amplitude of the 3D density fluctuations. As 
described in |3] the easiest way of obtaining the 3D power 
spectrum is to use the approximate relation between the 3D 
and 2D spectra according to eq. (|lHl . which for model power 
spectra works well at wavenumbers of order 10~^ kpc~^ and 
larger. It is also convenient to convert the final power spec- 
trum to the amplitude of the density perturbations using 
ea. (|12|l . The final results of this analysis are shown in Fig |14l 
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Figure 12. Top-left: raw XMM-Newton image; Top-Right: image with 15" circles around bright point sources removed; Bottom-Loft: 
surface brightness for a /3-model with /3 = 0.6 and core radius Tc = 245 kpc; Bottom-Right: XMM-Ncwton image divided by the /3-modcl. 
All images are 24' X 24' (or 648x648 kpc). Images are slightly smoothed with a Gaussian filter. 



This 3D amplitude characterizes typical variations of the 

density relative to the unperturbed value (i.e. — = Asd) 

P 

as a function of the spatial scale. The amplitude varies from 
~10% at scales of a few 100 kpc down to ~5% at scale of 
a few tens of kpc. At even smaller scales, pure statistical 
errors due to Poisson noise increase dramatically precluding 
tight constraints on the amplitude of the fluctuations. 



5 DISCUSSION 

The main conclusion one can draw from Fig: ll4l is that there 
are significant variations of the density fluctuations of the 
ICM on spatial scales ranging from ~500 kpc down to ~30 
kpc at the level of 5 to 10%. These fluctuations are rela- 
tively small over the entire range of wavenumbers where the 
existing data have sufficient statistics to make accurate mea- 
surements. Large-scale fluctuations are directly visible in the 
images divided by a suitable smooth model as excesses or 
dips in the surface brightness (Fig |12|) . although their ap- 
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Figure 13. 2D amplitude of X-ray brightness fluctuations (rela- 
tive to the underlying /3-model) as a function of the wavcnumber. 
Black - Chandra, red - XMM-Newton. The uncertainties shown 
are due to Poisson noise only and do not account for the stochas- 
tic nature of density fluctuations. The PSF correction (see i|4.4l l 
is taken into account. 

pearance does vary depending on the type of the smooth 
model. Small-scale perturbations are in the regime domi- 
nated by the Poisson noise and they are revealed only in a 
statistical sense and not individually visible in the image. 

There are several plausible physical mechanisms which 
can lead to the observed emissivity/density fluctuations. 
Among them are the following possibilities: 

• too simplistic model of the unperturbed gas density dis- 
tribution; 

• perturbations of the gravitational potential on top of a 
smooth global potential of the cluster; 

• entropy fluctuations caused by infalling low entropy gas 
or by advection of gas from one radius to another by gas 
motions; 

• variations of gas density and/or pressure associated 
with gas motion and sound waves; 

• metallicity variations; 

• presence of non-thermal and spatially variable compo- 
nents (bubbles of relativistic plasma and magnetic fields). 

We now briefly discuss these possibilities. 



5.1 Simplistic smooth model 

One of the principal uncertainties in calculating deviations 
of the surface brightness or density fluctuations lies in speci- 
fying the smooth underlying model. Our experiments with a 



k, kpc~' 

0.001 0.01 0.1 1 
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Figure 14. 3D amplitude of density fluctuations (relative to the 
underlying /3-modol) as a function of the wavenumber. Black - 
Chandra, red - XMM-Newton. The uncertainties shown are due 
to Poisson noise only and do not account for the stochastic nature 
of density fluctuations. The effect of the PSF has been taken into 
account in this plot. 



two-dimensional /3-model have shown that varying the cen- 
ter of the model, the core radii in two directions and /? 
within reasonable limits one can change the amplitude of 
the longest perturbations (400-500 kpc) by a factor of two, 
compared to a symmetric /3-model (see Fig[T5}. To illus- 
trate this statement three models were used in Fig llSI (i) 
spherically symmetric /3-model with /3 = 0.6 and Vc = 9'; 
(ii) two-dimensional /3-model with /3 — 0.62 and rc,x = 9' 
and rc,y = 8' and (iii) another two-dimensional model 
with /3 = 1.03 and rc.x = 15.9' and rc,y = 13.5' (from 
IVikhlinin. Forman. fc Jones! 1 19941 '). The variations of the 
amplitude at large scales is not surprising, since with more 
complex models we are able to attribute some of the struc- 
tures to the underlying model, rather than to perturba- 
tions. The amplitude of the perturbations at scales smaller 
than ~300 kpc remains unchanged. In spite of the model- 
dependent variations of the amplitude at large scales, the 
structures seen in the images divided by the /3-model are 
complicated enough, so it is probably fair to state that fur- 
ther reducing the amplitude of these structures would re- 
quire substantially more complex models. We concluded that 
a conservative systematic uncertainty - a factor of order of 
2 - should be attributed to large-scale perturbations, while 
at scales smaller than 300 kpc the systematic uncertainty is 
smaller. 
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Figure 15. 3D amplitude of density fluctuations relative to dif- 
ferent underlying /3-models (based on XMM-Newton data). The 
red points (the same as shown in Fie ll4l l correspond to a spher- 
ically symmetric /3-model with /3 = 0.6 and r"c = 9'. The blue 
points correspond to a two-dimensional /3-model with /3 = 0.62 
and rc,x = 9' and rc,y = 8'. The black points are for a model 
with /3 = 1.03 and Tc.x = 15.9' and rc,y = 13.5'. Depending on 
the parameters of the model the amplitude changes by a factor 
of 2 at scales above 500 kpc. At scales smaller than 300 kpc the 
power spectrum is virtually unchanged. 

5.2 Perturbations of the gravitational potential 

The gravitational potential of the cluster is affected by the 
presence of subhalos. In particular, in Coma there are two 
large cD galaxies - NGC 4874 and NGC 4889, with the 
X-ray surface brightness centroid lying in between these 
two objects. One can identify enhancements of the surface 
brightness at the positions of NG C 4874 and NGC 4889 
(|Vikhlinin. Forman. fc Jonej|l99^ ). although the exact am- 
plitude of these enhancements depends on the assumptions 
of the underlying cluster modeJf|. 

We can write the gravitational potential of the cluster 

as: 

V3 = V=o + SiA<^i, (23) 

where tpo is the smooth potential and A</?i is the perturba- 
tion, associated with the subhalo number i. To simplify the 
estimates we assume that for each of the two bright galaxies 
their potential can be written in an isothermal form over the 

® Note that when calculating power spectra the 15" circles 
around NGC 4874 and NGC 4889 have been excluded. There- 
fore the contributions of an AGN in NGC 4889 and cool coronae 
jVikhlinin. Forman. fc JonejlTggi ) in both galaxies are excluded. 



Figure 16. 2D relative amplitude of the surface brightness fluc- 
tuations (black line) caused by the presence of two components 
(representing gas concentration around two massive cD galaxies 
in Coma) on top of the global cluster emission. Parameters of 
the / 3-models for each galaxy and for th e main cluster are taken 
from lVikhlinin. Forman. fc Jone^ l|l994l ) . Note that the decompo- 
sition of the surface brightness into the three-component /3-model 
is not unique. The same set of exclusion regions (around compact 
sources, including nuclei of NGC 4874 and NGC 4889) as for real 
data was used for the simulated dataset. For comparison, the 2D 
amplitude measured from XMM-Newton data is shown with red 
crosses. 

interesting range of radii: 

Av3 = 1-^ In r-fC, (24) 

where Vc is the circular speed and r is the galactocentric ra- 
dius, C is an arbitrary additive constant. Further assuming 
that the gas is isothermal and ignoring possible effects of 
the galaxy motion with respect to the ICM, the relative gas 
density perturbation, induced by is took out the first 
minus sign in the equation but perhaps you prefer the other 
way 

J^^e^-l^(-) -1, (25) 

where fi ~ 0.61 is the mean particle atomic weight, rup 
is the proton mass, T is the ICM temperature, k is the 
Boltzmann constant and the overall normalization is ex- 
pressed in terms of the scaling radius, and the nor- 
malization constant C is c hosen so that Aip(rs) — 0. For 
NGC 4874 and NGC 4889, ICoccato. Gerhard, fc Arnaboldil 
(|2010l ) measured the line-of-sight velocity dispersion a of 
~ 283 km s~^ and ~ 266 km s~^ respectively at a distance 
of 50-60 kpc from their centers. For a spherical galaxy with 
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isotropic stellar orbits a reasonable estimate of the circular 
speed is then Vc ~ \/^a ~ 460 — 490 km s^^. Conservatively 
assuming Vc = 450 km/s and the ICM temperature kT — 7 
keV, the resulting radial dependence of the density pertur- 

n / \ -0.18 

5p 



bation is — oc — — 1. Further assuming that the 

P VsJ 

perturbation of the potential extends up to =150 kpc and 
is zero at larger distances, the RMS of the density variations 
induced by a single perturbation of the potential of this type 
can be estimated as follows: 



4nr dr 



1/2 



47rr^dr 



0.03, 



(26) 



where r^ax ~ 325 kpc is the characteristic size of the region 
which we analyze. To simplify this estimate we have assumed 
that the perturbation is located at the cluster center. This 
is a crude estimate, but it suggests that adding a single 
perturbation in a form of an isothermal halo with the above 
parameters is able to induce a volume averaged RMS of 
order 3%. 

Another possibility to make this estimate is to fit the 
observed enhancements in the surface brightness around po- 
sitions of the two bright galaxies and use the best fitting pa- 
rameters to estimate the RMS. To this end we used a decom- 
posit ion of the X-ray surface brightness into three /3-models 
from lvikhlinin. FormanT fc Jonej (|1994| ) - one for the main 
cluster and one for each of the two cD galaxies. Repeating 
the analysis for the three /3-models image divided by the 
/3-model for the cluster, we obtained a scale dependent rela- 
tive amplitude of the surface brightness fiuctuations induced 
by these galaxies (see Fig[T6]). For comparison we show the 
density fluctuations derived from the XMM-Newton image. 
Clearly, the amplitudes in the range of scales of order few 
100 kpc are in the ballpark of values predicted by this sim- 
ple exercise. This suggests that perturbations of the poten- 
tial can be responsible for a significant part of the observed 
perturbations. 



5.3 Turbulence 

Turbulence is another possi ble source of pressur e/density 
fiuctuations in the ICM fe.g. lSchuecker et al.|[20o3 ). Assum- 
ing an eddy with characteristic velocity v^ddy and adiabatic 
behavior of the gas one can relate the velocity and density 
perturbations as 



5p 
P 



2 7 itr 



(27) 



where 7 = 5/3 is the gas adiabatic index, ^ ~ 0.61 is the 
particle mean atomic weight, rup is the proton mass and 
M is the Mach number. For kT = 7 keV, the value of 
Vnddy ~ 450 km s~^ will induce 5.5% variations in the ICM 
density. Approximately (up to a factor of order unity) this 
corresponds to the ratio of the kinetic and thermal energies 
of the ICM. According to nu merical simulations this is a 
perfectly reasonable value (e.g . Sunvaev. Norman, fc BrvanI 



if numerical simulations can accurately predict the charac- 
teristics of the turbulent eddies (size and velocities), since 
the Reynolds numbers achieved so far are not very high (and 
we in fact do not know the effective viscosity of the ICM). 

The assumption that a significant part of the observed 
fiuctuations is due to perturbations of the gravitational po- 
tential (see i]5.2p would imply that the level of turbulence is 
modest. For instance one can use Fig |16l to estimate the am- 
plitude of fiuctuations on top of the contribution induced by 
gas concentration around the two massive cD galaxies. These 
excess fiuctuations of density have an amplitude of order 2- 
4% on scales ~30-300 kpc, which corresponds to a velocity 
amplitude of ~ 300 — 400 km s^^ according to eg 1271 This in 
turn implies that the contribution of turbulent pressure sup- 
port is similarly low. Accounting for the gas entropy fiuctua- 
tions considered in section i]5.5l mav lower this estimate even 
further. Of course accurate estimates of the turbulent pres- 
sure support by this method are difficult. Future observa- 
tions w ith X-ray micro-calorim eters, initially from ASTRO- 
H fe.g. iTakahashi et al.llioiol '). should be able to constrain 
both the velocities and sizes of the energy-containing eddies 
(e.g. IZhuravleva et al.llioill '). For now, turbulence remains a 
plausible contributor to the observed surfac e brightness fluc- 
tuati ons, as well as pressure fiuctuations (jSchuecker et al.l 
l2004h in Coma. 



5.4 Sound waves generated by the turbulence 

It is interesting to estimate a possible contribution of sound 
waves to the observed surface brightness variations. The ad- 
vantage of sound waves is that they represent a linear rela- 
tion (rather than quadratic) between the gas velocity and 
the pressure/density perturbation: 



5p 
P 



(28) 



where is the sound speed. E.g. to produce 5% variations 
in density one needs u ~ 70 km s~^. Sound waves cross the 
region of interest 7? ~ 350 kpc over ts = R/cs- If one as- 
sumes that the sound waves are generated by the turbu- 
lence, then the rate of so und-wave gene r ation per unit vol- 
ume (ergs s^^ cm^'^) is (|Lighthill|[l95^ : iLandau fc Lifshit j 
§75) 



"sddy 



(29) 



where I is the characteristic eddy size. The typical amplitude 
of the gas velocity due to sound waves in the region of size R 
(assuming that turbulence generates sound waves uniformly 
over the region with volume ~ i?^) can be estimated as 



/ ets = V, 



eddy 



' y ~ 30 km s ^ 



(30) 



l2003l : lLau. Kravtsov. fc Nagaill2009l ) . It is not clear however 



where we assumed Veddy = 450 km s~^ and / — 100 kpc. 
Thus sound waves generated by turbulence, are not very 
effective in producing density perturbations unless the tur- 
bulent velocities generating the sound waves are larger than 
450 km s~\ Of course the above argument provides only an 
order of magnitude estimate. 
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5.5 Entropy variations 

Another natural mechanism, which may cause density vari- 
ations is the variation of entropy/temperature. We distin- 
guish them from potential perturbations, turbulent pulsa- 
tions or sound waves, since entropy variations can exist with- 
out associated strong pressure perturbations (to lowest order 
in Mach number). We assume that in these perturbations 



Sp 
P 



(31) 



where T is the gas temperature. 

One can identify two different modes of entropy pertur- 
bations. The first mode is associated with the low entropy 
gas, infalling into the system. One of the features, seen in the 
Coma images, namely a c hain of galaxies apparently infalling 
into Coma from the East (|Vikhlinin. Forman. fc Joneslll997 ) 
is associated with colder gas, which is clearly seen in the 
hardness ratio map (e.g. [Neumann et al.ll2003l ). In Fig|12l 
this cold stream appears as an enhancement in the surface 
brightness relative to the cluster /3-model. The amplitude 
and the area covered by this feature are comparable with 
the perturbations, identified as a potential perturbation in 
8^5.21 We therefore can expect a similar impact of this feature 
on the final estimate of the density perturbations. 

Another possibility to create entropy fluctuations is the 
advection of gas with different entropies by gas motions in 
turbulent eddies. Assuming an unperturbed isothermal at- 
mosphere of the cluster with density following a /3-model, 
we can estimate the eddy size needed to create the required 
density contrast. When a gas lump is advected from the ra- 
dius r, where the gas pressure is P, to the radius r -\- 5r, 
where the pressure is P — 5P, and it expands adiabatically, 
the corresponding density change is 



Sp 
P 



15P_ 
7"P 



7 



-3/3 



6r 



1 + 



(32) 



Near the core radius r ~ sl 10% change in radius causes a 
~5% change of density. Thus eddies of size 5r ~ O.lVc ~ 25 — 
30 kpc are needed. As in the case of turbulent fluctuations, 
future measurements of eddy amplitudes and sizes would 
provide crucial constraints on the contribution of this mode. 



5.6 Metallicity variations 

In the 0.6-2.5 keV band the change of metalli city from 0.5 
solar to 1 solar (for kT = 8 keV APEC model, ISmith et al] 
(|200ll )') causes a 5.6% increase in the count rate (for XMM- 
Newton). The change of density, causing the same increase 
of the flux is 5.6/2=2.8%. Thus very strong variations of 
metallicity are needed if they are to make a significant con- 
tribution to variations of emissivity. 



5.7 Bubbles of relativistic plasma 

Observations of cool-core clusters revealed rich substructure 
caused by the pr esence of AGN-inflated bubbles of relativis- 
tic plasma (e.g. IChurazov et al.ll2000l ). In terms of X-ray 



surface brightness perturbations, the bubbles reveal them- 
selves as depressions in the surface brightness with the am- 
plitude set by the fraction of the line-of-sight volume oc- 
cupied by the bubbles. In cool-core clusters, the fraction of 
the volume occupied by bubbles can be large - of order 10% 
or more. The RMS variations of the density perturbations 
will be high accordingly. In Coma there is no compelling 
evidence for widespread bubbles of relativistic plasma, al- 
though one of the cD galaxies - NGC4874 does contain an 
active ra dio source with a pair of bubbles on each side of the 
nucleus (|Sun et al.ll2005l l . Future low frequency radio obser- 
vations might reveal radio bright regions (associated with 
low energy electrons contained in bubbles). If these regions 
correlate with the X-ray surface brightness depressions, this 
would be a key indicator in favor of this mechanism. For now 
we can only state that in Coma, radio bubbles are a possi- 
ble source of X-ray surface brightness fluctuations, although 
there exists no compelling evidence that this mechanism is 
really important for Coma. 



6 CONCLUSIONS 

X-ray surface brightness fluctuations in the core (650 x 650 
kpc) region of the Coma cluster are analyzed using XMM- 
Newton and Chandra images and converted into scale- 
dependent amplitudes of volume-filling ICM density pertur- 
bations. The characteristic amplitudes relative to a smooth 
underlying density are modest: they vary from 7-10% on 
scales of 500 kpc down to ~5% at scales ~ 30 kpcQ. At 
smaller scales, projection effects and the PSF (for XMM- 
Newton) smear the density fluctuations by a large factor, 
precluding strong limits on the density fluctuations in 3D. 
Nevertheless making deep (~ Msec long) Chandra exposures 
can still be useful in pushing down the constraints for scales 
smaller than 30 kpc. 

Several physical effects can contribute to observed fluc- 
tuations of the surface brightness/density. The most likely 
contributors at large scales are the perturbations of the grav- 
itational potential by massive cD galaxies and entropy vari- 
ations due to infalling cold gas. We do see in raw images the 
substructure which can be plausibly associated with these 
mechanisms. 

Other effects can contribute to the observed fluctua- 
tions, including variations of density in turbulent eddies, ad- 
vection of the low/high entropy gas by turbulent eddies, gen- 
eration of sound waves by turbulence, bubbles of relativistic 
plasma and metallicity variations. Future high energy res- 
olution X-ray missions (e.g. ASTRO-H) and low frequency 
radio observations (e.g. LOFAR) could help differentiate be- 
tween these possibilities. At present we can consider the de- 
tected variations as conservative limits on the amplitude of 
the density /emissivity fluctuations caused by any of these 
effects. 

Despite a variety of possible origins for density fluctu- 
ations, the gas in the Coma cluster core is remarkably ho- 



^ for comparison the mean free path for Coulomb collisions in 
the core of Coma is of order 5 kpc 
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mogeneous on scales from ~ 500 to ~30 kpc, with density 
variations relative to a plausible smooth model smaller than 
10%. 
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